The numerical modeling of lava dome evolution at volcán de Colima using VOF and SPH methods

Lava flows from extrusive volcanic eruptions can have catastrophic consequences both for human life and the environment. Modeling such situations is an important scientific problem. The main driving forces in evolution of mentioned lava are gravitational forces, viscous friction on surface spill, processes crystallization molten rocks into plateau, tubes, domes. In this paper, mathematical model eruption includes Navier–Stokes equation, incompressibility phase transfer as well corresponding initial boundary conditions. Mathematical models considered compared within Euler (Volume Of Fluid – VOF) Lagrange (Smooth Particle Hydrodynamic SPH) formulations. ANSYS Fluent, OpenFOAM, SPlisHSPlasH packages were used computer simulation. Computer simulation algorithms problem implemented C++ language. Numerical modeling a real dome formed at Colima volcano (Mexico) February–March 2013 was carried out. For experiment, information about dynamics growth, collected during eruption, used. It shown how approach makes it possible to establish dependence morphology rheology highly fluid intensity outflow.